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Abstract 

We summarize recent attempts to extract characteristics of non-perturbative 
vacuum structure from lattice measurements of the gauge invariant field 
strength correlator. As an alternative to cooling, we propose to apply the 
renormalization group (RG) smoothing method in lattice studies. For pure 
SU(2) gauge theory we present magnetic and electric correlation lengths 
and condensates related to various correlators over a temperature range of 
0.7 T dec < T < 1.9 T dcc . 
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I. INTRODUCTION 



Quantum Chromodynamics (QCD) owes its confining property to non-perturbative struc- 
tures of its groundstate which are, after many years of research, not yet uniquely identified, 
(for a recent overview of competing pictures see e.g. In a more pragmatic attitude, 

the Stochastic Vacuum Model (SVM) 0] describes the vacuum structure in terms of a few 
vacuum correlation functions and points out which are related to confinement. That only 
few correlators should characterize the vacuum is the consequence of a truncation of the 
hierarchy of correlation functions, essentially due to a postulate made within the SVM that 
the gauge fields are Gaussian random processes in Euclidean space-time. In this framework, 
the correlators enter, and relate to each other, various phenomenological predictions. Par- 
ticularly interesting (exhausting, in the case of pure gluodynamics) is the gauge invariant 
gluon field strength correlator in Euclidean space (and its continuation to Minkowski space). 
According to the SVM, part of the Euclidean field strength correlator at zero temperature 
is related to the confinement property of the QCD ground state ||. For heavy quarkonium 
states, the full correlator describes the effect of the gluon condensate on the level splitting in 
detail (see Ref. || and references therein). Wick- rotated to Minkowski space, the correlator 
serves to parametrize high energy hadron-hadron scattering (see e.g. Ref. H). 

Regardless of how this truncation exactly works, these correlators incorporate the un- 
known non-perturbative mechanisms and can be considered as an implicit parametrization 
of the important gluonic vacuum field modes. If they were known, this information could be 
readily used to constrain the parameters of more explicit models of vacuum structure. The 
central role, both theoretically and phenomenologically, strongly suggests to extract the field 
strength correlator (and other vacuum correlators) from a first principle lattice simulation. 
It is interesting to notice that, although the SVM approach has been proposed more than a 
decade ago by Dosch and Simonov 0, on the lattice the field strength correlator has only 
been studied by the Pisa group. Apart from earlier results for SU (2) || and SU (3) pure 
gauge theory, there are now precise results available, from an analysis using cooling, over an 
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interval of distances 0.1 ... 1.0 fm for pure SU(3) gluodynamics @||, for QCD with dynam- 
ical (Kogut-Susskind) fermions |L0| (four flavors, for two values of mass) over an interval 0.4 
... 1.0 fm, and, over the same interval of distances, for pure SU(3) Yang-Mills theory very 
near the deconfming temperature M. 



Lattice data for T = obtained and fitted in Refs. P|^TU[ give strong hints in favor of 
dominance of locally selfdual or antiselfdual field strengths. These data have been analyzed 
TT| from the point of view of the random instanton liquid model |l^JT3[] . Here the following 



picture arises: a rather dilute liquid of instantons (with parameters which, in the case of full 
QCD, happen to be roughly consistent with parameters obtained from phenomenological 
analyses) can explain the field strength correlator, including the part which is believed to 
be responsible for confinement. The relevance of this statement for confinement is far from 
clear. Absolute confinement is, strictly speaking, a property of pure Yang-Mills theory alone. 
At this place let us mention that a much higher instanton density is extracted from the field 
strength correlator in the case of pure Yang-Mills theory than for full QCD ||1 1|| . 

On the other hand, according to a widespread opinion a random instanton liquid, how- 
ever dense it might be, does not lead to an asymptotic string tension. This is obvious for 
calculations of the direct effect of the (unperturbed) instanton liquid on the static Q-Q 
potential via the Wilson loop. While the instanton potential rises linearly at intermediate 
distance, this does not correspond to the known string tension |14| if the phenomenolog- 



ical density and size are assumed to hold. Let us note, however, that neither instanton 



suppression by external charges or fields nor instanton interactions (see e.g. [15|]) have been 
taken into account in these estimates. It has even been questioned |16| whether the available 
lattice evidence for confinement is conclusive at all. In short, in spite of the possibility of 
an instanton interpretation of the correlator, the role of instantons for confinement is still 
controversial. The answer depends on details of the modeling of the instanton liquid, in 
particular on the density and the shape of size distribution (see e.g. Refs. [|17|,|I8|1). 

There is an indirect route (not completely worked out yet) leading from the instanton 
liquid to confinement. It has been demonstrated that under certain circumstances an instan- 



ton liquid induces the condensation (percolation) of magnetic monopoles (see Ref. |1| and 



earlier references therein) which should confine via a dual Meissner effect. A corresponding 
effective theory of a dual Abelian gauge field, interacting minimally with the (condensed) 
monopole field (Dual Abelian Higgs Model JOJ) is characterized by the mass of the dual 
"photon" and a finite Higgs correlation length. Both quantities together determine PUj the 



gauge invariant non-Abelian field strength correlator which we are interested in. Therefore, 
in a more fundamental way, the two structure functions of the correlator at T = are 
inherent to an infrared effective dual formulation of QCD (see also f21~H ). 



Similarly to T = the field strength correlator at finite temperature encodes information 
on generic non-perturbative structures of non-Abelian gauge field histories dominating the 
QCD partition function. In a similar spirit as of Ref. the correlators obtained near 
the deconfinement transition in pure SU(3) gauge theory |HJ have been analyzed from the 
caloron (finite T instanton) point of view p2| . Such a description turns out to be possible 



only below the deconfinement temperature Td ec . 

As any measurement involving the field strength, such lattice calculations are far from 
trivial. The cooling method has been used in most of the more recent studies. To our 
opinion, the implications of this method for the intended measurements have not been 
discussed thoroughly enough so far. In this paper, we are going to apply a different technique 
to improve the signal for extracting the field strength correlator from lattice gauge field 
ensembles. Our method is based on a (classically) perfect action and uses a renormalization 
group motivated smoothing procedure |2^j2^]. We shall give a presentation of first results 
obtained in this manner for finite temperature SU(2) pure gauge theory. In contrast to 
T = 0, a similar splitting of correlation lengths related to the different structure functions is 
obtained with both techniques. Some possible differences should be closer examined whether 
they are due to the other gauge group SU{2) used here, compared to SU (3) analyzed with 
cooling. If not then, cooling would have a systematic effect concerning the relative strength 
of the structure functions D and T>% at high temperature. 

The paper is organized as follows. In Sec. II we give an introduction to the concept of 



field strength correlators for zero and finite temperature; the correlators and their suitable 
decompositions into the structure functions are defined there. In Sec. Ill we briefly discuss 
the potential problems of the cooling method and recall in short the fixed point action 
renormalization group (RG) smoothing technique. In Sec. IV we describe our results, and 
conclude with a discussion in Sec. V. 

II. FIELD STRENGTH CORRELATOR 

The two-point field strength correlator is a non-local generalization of the gluon conden- 
sate, involving the field strength tensor at separate points : 



where G w = T a G a w is the field strength, with T a being the generators of the algebra of 
the color SU(N C ) group in the fundamental representation. S(xi,X2) is the (path ordered) 
Schwinger line (parallel transporter) 



with = T a A^. It is necessary to include the non-Abelian Schwinger line phase factor in 
order to ensure gauge invariance. Usually the straight-line path from x\ to x 2 is chosen to 
evaluate S(x 1 ,x 2 ). The correlator describes the coherence of the gauge field (modulo gauge 
transformations) over some coherence length. This notion replaces the idea of "domains" 
originally imagined to form the gluon condensate. Completely contracted, it is the local 
gluon condensate defined as 



= 3 2 (0|Tr {G w (0)S(0,x)G w (x)S l (0 ) x)} |0> , 



(1) 




(2) 



(oe a 




(3) 



which is obtained in the limit \xi — x 2 



from the regulated correlator. 



A. Zero Temperature 



At T = 0, Euclidean 0(4) invariance admits the decomposition (for x = X\ — x 2 ) 
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Vp, pu(T {x) = (8 ^8 - Sfj^Spv) [V(x) + Vi{x)\ 

, <9£>i(x) 

-|- [X pXyO p a X pX a O py ~\- X pX a O p V X pXyO p a ) Qrgl ' \ I 

where T> and T>\ are two invariant functions (of Euclidean modulus \x\) containing informa- 
tion on the vacuum. The tensor structure related to T>\ does not contribute to confinement, 
while T> 7^ is possible only in non-Abelian theories or in Abelian theories with monopoles. 

It is convenient, in particular on the lattice, to define two other vacuum structure func- 
tions T>\\(x) and a D±(x), according to whether the plaquettes defining the field strength 
component in question have one direction parallel to the 4-vector between the two points of 
observation. These functions are composed of T> and T>\ as follows 

V\\ = V + V 1 + x 2 —^ , 
ox z 

V ± = V + V 1 . (5) 

The two basic vacuum structure functions at T = have a non-perturbative and a pertur- 
bative part according to the notation 



V(x) = V np (x) + ^e 



^-N/Ao 

x\ 



V 1 (x)=V7(x) + ^e-^ , (6) 

fI 

Since the first structure function is closely related to confinement, the non-perturbative parts 
of the two functions are suggestively written as 

V«*{x) = ^G 2 V(x) , 
o 

V7(x) = ^(1-k)G 2 V 1 (x) , (7) 
o 

where G 2 has to be identified with the gluon condensate (||) and 

V(0)=V 1 (0) = 1. (8) 

The so-called non-Abelianicity k = T> np (0)/ (T> np (0) +P™ P (0)) represents a certain fraction 
of the gluon condensate related to confinement. So far the functions T>(x) and T>\{x) mostly 
have been fitted as exponential functions [|?HT0|] of the modulus \x\ 
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f)(x) = e" |x|/Ao and =e" |a;|/Al , (9) 



although more physically motivated forms are available which fit the data as well [2S|. 

In the SVM the string tension can be approximately related to the correlation length Ao 
(coherence length of V) and the gluon condensate @, 

o = - J d 2 xV(x) 

(10) 



3 3 

— kG 2 [ dx 2 V(x) = — k G 9 M . 
12 J K ' 6 



In the perturbative regime (at short distance) both T> and T>\ behave as l/|x| 4 . One- 
gluon exchange contributes only to the latter [2B| giving T>\ = g x ^ ■ One- loop and higher 
order perturbative contributions to the string tension integral flTo|) are cancelled by higher 



correlator contributions 27 



Purely selfdual or antiselfdual non-perturbative background fields contribute only to T> 



(see Ref. [0,|T|]). In this case, it is the overlap of different selfdual or antiselfdual configu- 
rations (higher order contributions in an expansion in the density |TT|) which generates the 
non-perturbative part of T>\. The observation that T> » first made for T = lattices, 
supports the idea that strong fields are dominantly selfdual or antiselfdual. 

It was only for T = that the newer data have been analyzed by the authors themselves 



As short distances became observable at higher (5 it has been confirmed that the 
functions T> and T>\ contain a perturbative part which could be subtracted by a fit like (^|). 
For the case of quenched SU(3) gauge theory the two correlation lengths have been found 
equal, Ao ~ Ai ~ 0.22 fm. The data for full QCD show that the correlation length Ao 
becomes smaller with increasing quark mass towards the quenched case. For the chiral limit 
an upper limit for the correlation length can be estimated as Ao ~ Ai ~ 0.45 fm. The gluon 
condensate G2 as extracted from the correlator according to (HI?]) is much larger in pure 
Yang-Mills theory than any phenomenological value. 



1 For quenched SU(2) gauge theory, a first analysis || without cooling has found Aq ~ 0.17 fm. 
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In Ref. jn| the field strength correlator data at intermediate distance have been used 
to extract the two parameters characterizing the random instanton liquid model [[12] which 
gives a good description of most of the hadronic correlators at zero temperature. Q The 
size distribution of this model is approximated as d(p) = 5(p — p), so density n 4 and size p 
are the only parameters. The investigation in Ref. |TT| was carefully dealing with the path 
dependence of the Schwinger factor in the instanton background. For a straight line path, a 
reasonable description over the measured range of distances has been achieved with param- 
eters roughly compatible with phenomenological estimates where this is possible: density 
n.4 = 0.5 fm -4 and size p = 0.44 fm for full QCD (with four flavors), to compare with 
n 4 = 4 fm -4 and p = 0.3 fm in pure SU(3) Yang-Mills theory. Second order terms in n 4 give 
subleading contributions to T> and are the leading ones in T>\. As long as deviations from 
selfduality or antiselfduality are entirely due to II, II and II overlap, this analysis has led 
to the conclusion that T>\ < is characteristic for the instanton liquid at T = 0. Although 
this was not the case in the original fit given by the authors, it has been shown that the 
quenched data are compatible with a small negative T>\ as well. With the second order 
terms taken into account for pure Yang-Mills theory, the best fit leads to a change towards 
smaller n 4 = 3.4 fm -4 and p = 0.27 fm. The smallness of the correction a posteriori justifies 
the dilute gas approach to pure Yang-Mills theory and even more so for full QCD. As a rule 
one can summarize that instanton radii estimated in this way are 20 % to 30 % bigger than 
the (exponential) correlation length of T>. 



2 Strictly speaking, with the classical shape of an instanton in the trivial vacuum the asymptotic 
behavior at x 2 — > 00 of the correlator T>(x) can not be described. The inclusion of the interaction 
with the vacuum medium into the "classical" Euclidean field equation results in a modified instan- 
ton profile far from the center [29]. For a recent interesting attempt to find a modified instanton 
solution and to use it for the correlator see Ref |3C]. 



8 



B. Finite Temperature 



At finite temperature, the 0(4) space-time symmetry is broken to the spatial 0(3) sym- 



metry and the number of independent correlators is enlarged [31|. All of them separately 
depend on |xi — x 2 | and (x 4 )i — (x 4 ) 2 . We consider 



(0\Tr{E i (x 1 )S{x 1 ,x 2 )E j {x 2 )S\x 1 ,x 2 )}\0) 



r f^F ^F 29T>f\ 9Vf 

5 ij (v E +V E + x 2 —^)+x i x j 



(11) 



Y dx'j J ' " J% 3 dx 2 

where Ei = Ga is the electric field, and 



(0\Tr{B i (x 1 )S{x 1 ,x 2 )B j (x 2 )S^{x 1 ,x 2 )}\0) 



ki [V B + V B +x 2 — L )-x i x j 



(12) 



^ dx 2 J ~ % 1 dx 2 

where Bi = ~ eijkGjk is the magnetic field. Let us assume that x = (0,0, |a?|,0) and define 
longitudinal and transversal correlators for electric and magnetic fields 

Vf = V E + V? + x 2 ^ (i — j — 3 in eq.(ll)) 



VI = V E + Vf (i = j = 1, 2 in eq.flnD) 

(13) 

Vf = V B + V B (i = j = 3 in eq. ([12])) 
f)V B 

V B = V B + V B + x 2 ^ (i=j = l,2 in eq.©). 

In addition there is an electric-magnetic correlator which will not be considered here. 

As a first step, in Ref. |22| an exponential plus perturbative fit in terms of V E,B and V E ' B 



has been tried to describe the finite T quenched SU(3) correlator measured in || with cooling 
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(which is known at distances > 0.4 fm). For the magnetic correlators, immediately below 
the deconfining transition, at a temperature T = 0.978 T^ ec , the two structure functions 
have clearly different correlation lengths, = 0.19 fm and Af = 0.53 fm. The first one 
slowly increases across the transition, while the second goes through a maximum there. For 
the two electric structure functions at the temperature T = 0.978 T^ ec , different correlation 
lengths A^ = 0.17 fm and Af = 0.41 fm have been found, too. These fits show a tendency 
for both electric correlation lengths to drop suddenly at the transition. Unfortunately, 
there is no cooling data for pure SU (2) at finite temperature that would allow a more 
direct comparison. A description using the Harrington-Shepard caloron solution, i.e., the 
periodic instanton |32|,|33[], is the main objective in Ref. p2[ . Below the deconfinement 
temperature, where T> B m T> E and T)f ~ T>f has been found and where T> EB » Z>f ' B , a 
satisfactory description of magnetic and electric correlators was possible using these selfdual 
or antiselfdual configurations while the small differences could be attributed to overlapping 
two-caloron (or caloron-anticaloron) configurations. The ratio T>f /T> E strongly increases 
above the transition. 

After describing our smoothing method in the next section, we will discuss the correlation 
lengths for the different finite temperature structure functions according to exponential 
Ansatze similar to what has been done in Ref. |E| for the cooled SU(3) data, making 
combined fits to electric and magnetic data separately, in terms of T> E and T>f and in terms 
of T> B and T> B , respectively. We will also discuss a fit which directly applies to the raw data 
for Pi| and T> E simultaneously, as well as for T> B and T> B . In this case we will separate the 
data into D and D\ afterwards, using the obtained form of the fit. It should be emphasized 
already here that, due to our smoothing method, these functions are obtained essentially 
free of a perturbative contribution. 



III. LATTICE METHODS 
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A. Smoothing vs. Cooling 



The Pisa group [0-10] has used the cooling method in order to measure the field strength 
correlators at T = 0, both without external sources (in vacuum) and in the presence of an 
external Wilson loop [34], [] and at temperature near T^ ec . For the measurement, similar to 



the search for instantons, it turned out necessary to eliminate ultraviolet fluctuations. It has 
been argued that, due to the diffusive nature of cooling, the correlator to be measured at a 
distance da (a = lattice spacing) would not be affected before a number n coo i oc d 2 of cooling 
iterations would have been performed. Practically, one measures the expectation values 
defining TD\\{dd) and T>±(da) as a function of n coo i- The actual correlator is then defined as 
the plateau value vs. n coo i- Although it has not been presented in detail, experience seems 
to justify the expectation that this plateau is reached safely before \x\ — x 2 \ happens to 
fall within the diffusion radius d oc y/n coo \. In another context it has been shown that this 
amount of cooling amounts to a renormalization of coupling, which is expected to show up 
in the perturbative signal contributing to the correlator. 

If cooling is technically necessary to improve the signal-to- noise ratio, one could be 
concerned to what extent it distorts the non-perturbative modes that existed in the true 
vacuum, i.e., to what extent the correlation lengths and the strength of the correlators are 
affected by the procedure itself. Since the (exponential) correlation length of the correlator 
T> is even smaller than the phenomenological instanton size it is useful to recall what effect 
cooling has in the attempt to fix instanton size and density from analyzing the topological 
density of a configuration. For this purpose, the cooling technique nowadays is superseded by 
adopting two, essentially opposite extreme ways to control that effect. The first prescription 
requires to observe the average size as function of n. coo i, extrapolating it backward to zero 



cooling steps |3Bj in order to obtain the would-be instanton size in vacuum. The other uses a 



3t 



Bali et al. [35 have analyzed the correlator in the presence of external charges, analyzing 



quenched lattice data for heavy quark forces adopting certain factorization assumptions from SVM. 
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global stopping criterion [|37| requiring cooling to be stopped when the density of instantons 
becomes equal to the (much better known!) topological susceptibility. This second method 
implicitly adopts a non-interacting instanton gas picture, which differs principally from the 
instanton liquid. It is recognized by now that one has to control cooling since instantons, once 
they become discernible by cooling or smearing, are steadily growing, becoming more and 
more classical, while small size instanton-antiinstanton pairs annihilate. Improved cooling 
]38| aims to stabilize instantons with p > p m i n , still suffering from the latter problem. 
Presently, there is no method agreed upon, such that the instanton size cannot be uniquely 
defined on the lattice |39|. Each of the methods seems to have a smoothing radius of the 
lattice fields (depending on n cool ) which should be well below the distances at which one 
measures the correlator. Q 

During the last two years a method has been advocated p3| , |24|| that avoids propagating, 
iterative procedures to suppress UV fluctuations. Alternatively, it uses first blocking and 
subsequent inverse blocking with respect to a perfect action in order to accomplish the nec- 
essary smoothing of equilibrium configurations, while it preserves their long range structure 
beyond the smoothing scale. So far, this method has been used exclusively to analyze the 
instanton |23|] and monopole structure PB| , |2^J^T[] of Monte Carlo generated configurations. 



As for monopoles, smoothing wipes out small monopole loops which are irrelevant for con- 
finement (similar to blockspin renormalizations directly applied to the magnetic currents 

EH)- 



B. Smoothing with a Perfect Action 

To carry out a theoretically and technically consistent analysis, one has to implement a 
fixed point (FP) action. The simplified fixed-point action |^3[ for SU (2) has been improved 



4 Very recently a local stopping criterion for cooling has been proposed [[KJ that seems to accom- 
plish a tunable, well-defined smoothing radius. 
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recently, and a suitably truncated version [FT] lias been obtained. It is also parametrized in 



terms of only two types of Wilson loops, plaquettes (type C\) 

Ud = U Xtft7 u = U X7l j,U x+IM)V U x+VjlJ U^ v (14) 
and tilted 3-dimensional 6- link loops (type C 2 ) of the form 

Uc 2 = U X)fl ^ x = U x ^U x+ ^ v U x+ ^p i xU x+v+ ^JJ^ +XtV U X)X , (15) 

and contains several powers of the linear action terms corresponding to each loop (of both 
types) that can be drawn on the lattice 

Sfp(U) =EEE4,J)(1 - ~ Tr U Ci ) j . (16) 

i=l C t j=l / 

The parameters of this action as used in this work are reproduced in Table I. A comparison 



reported in detail in |24| shows that the new parametrization is near to perfect. 

With this action at hand one can now proceed to perform renormalization group trans- 
formations on the lattice. Suppose we have a fine lattice of links U and a coarse lattice of 
links V covering the same physical volume. Being a classically perfect action, Spp can be 
evaluated on both lattices and must satisfy the following condition 

Minu (S FP (U) + k T(U, V)) = S FP (V) , (17) 

configuration by configuration over a representative ensemble of equilibrium gauge field 
configurations V, each in one-to-one correspondence to a fine configuration U saturating 
the lower bound on Sfp(U) + k T(U,V) provided by the right hand side. T(U,V) is a 
certain non-negative functional related to the blockspin transformation (see |43| , p3| ). Eq. 
(pi) can be read as the saddle point equation for an integral over fine configurations defining 
the blocked action Spp(V). Neglecting the one-loop corrections is justified in the classical 
limit for (3 —>■ oo. The parameter « = 12.0 (not to be confused with the non-Abelianicity 
factor mentioned above) has been fixed in the context of the adopted blocking scheme by 
an optimization of the locality properties of the quadratic part of the action. 
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Blocking means a mapping from the link fields U on the fine lattice to link variables V 
which describe the same configuration on the next coarser length scale. This mapping is 
accomplished by the following construction: 

T> _ Mock jj jj 

+ C 2 \J^x,uU x +i>,iiUx+i>+(i,^U^.2fi,v U x _0^Ux-v,ij,U x -P+{i^U x -j)+2(i,i 



The coarse link variable V x ^ is obtained by normalization to a SU(2) group element V 3 



X,fl 



V x JJdet (V x ^). The blocking parameters c b 2 lock = 0.12 and cf ocA = 1-6 c h 2 lock have been 



optimized earlier in Ref. |^5[ under rather general conditions. 

Inverse blocking is a mapping V -> C/ SM where the new link variables U = U are 
implicitly defined as the set of links providing, for the configuration encoded in V, the 
"smoothest interpolation" on the next finer grained lattice. Practically, the link field JJ SM 



is found by minimizing the extended action on the left hand side of Eq. (TT). This is not 
a truly classical field on the lattice but it is smoother than any quantum field U that could 
give rise to the coarse level field V. Technically speaking, k plays here a role similar to a 
Lagrange multiplier defining the strength of the constraint imposed by the coarse grained 
configuration. 

The SU (2) pure gauge field configurations to be analyzed below have been generated on 
a 12 3 x 4 lattice at various f3 with the (classical) perfect action. In Ref. |24] we have localized 



the deconfinement phase transition at j3 c = 1.545(10). Corresponding to the second order 
character of that phase transition, the intersection of Binder cumulants of Polyakov loops 
for different 3d lattice volumes had been used for this purpose. Although the simulation was 
done on the fine lattice, the Binder cumulants built from V lattice Polyakov loops give the 
same f3 c . While the Polyakov line correlators become critical at the transition, the various 
field strength correlation lengths remain finite across the transition. 
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C. Measuring the Field Strength Correlators 



Here we demonstrate how we can get immediate access to the non-perturbative part of 
the field strength correlators applying the RG smoothing method. With some reservations, 
the results can be compared with cooling analyses (unfortunately available only for 577(3)). 
The field strength components g a 2 G fll/ (x) are defined as the average of the antihermitean 
part ^ (U x ,fi,u — U^u) °f the f° ur equal-oriented "clover leaf" plaquettes in the [iv plane 
which are open at the lattice site x: 

g a 2 G^(x) = ——7 {U XtlttV + V x ,v-n + U x ,-n-u + U x -v,n ~ h.c.) . (18) 

o I 

For the purpose of the following study we take advantage of the fact that in the result of 
smoothing, similar to cooling, the renormalization constant becomes Z ~ 1 for the gluon 
fields. 

We have measured only equal-time correlators, i.e. with a timelike argument (x 4 )i — 
(x 4 )2 = 0. For lattice sites x\,x 2 separated not along one of the spacelike axes, the choice 
of the Schwinger line is not unique. In that case averaging of expression ([]]) has been 
performed over a random selection of ten shortest paths connecting x\ and X2, each path 
defining an independent transporter S{x\ 1 X2). Naturally, for Pjf, T>¥, Df_, and T>f_, the 
field components Eh, Bn, E±_ and parallel or perpendicular to the 3- vector X\ — x 2 have 
been appropriately projected out of the field strength tensor G^. 

IV. RESULTS 

The raw data for (3 = 1.4 and f3 = 1.5 (confinement), for f3 = 1.55 (just above decon- 
finement) and three temperatures in the deconfined phase (fl = 1.6, (3 = 1.7 and (3 = 1.8) 
are shown in Fig. 1, together with an ad hoc fit to be discussed below. The two correla- 
tion functions T>f_ and T>^ have to be understood as averages over two correlators involving 
mutually parallel components of E or B orthogonal to X\ — x 2 . For better visibility the 
sum is shown in Fig. 1. The data shown are averaged over correlators with exactly equal 
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distance. Even then the data still would show an eventual rotational non-invariance (in the 
three-dimensional distance) if rotational symmetry would be seriously broken. 



A. Two- Exponential Fit 

At first, analogously to the procedure of DiGiacomo et al. applied to their T = cooled 
data |8|-|T0|, we have adopted the decomposition (|13|) into the basic structure functions for 



the E and B correlators separately, here however without a perturbative part which is 
suppressed by the steps of blocking and inverse blocking (see Sec. IV.B). We have assumed 
the following form for the basic E and B correlators 

V{x) = A)e HX|/Ao , 

(19) 

Pi(aj) = A,e- m/Al . 

Only correlator data at distances d > 2 a have been included into the fit, in accordance to 
our understanding that fields below that scale are interpolated in a smooth way and taking 
the extendedness of the field strength operator fll8[) into account. The four corresponding 
parameters are given in Table II to allow a comparison with the amplitudes and exponential 
correlation lengths obtained for T « T& cc by the cooling method for SU (3) [[J as analyzed 



in Ref. 22 



First, we mention that all structure functions are found positive at these temperatures 
in both cases. This property is corroborated also in the second way of fitting to be discussed 
below. Generally, at all temperatures, the two-exponential fit does not properly describe 
the decay of T> B as measured here by the smoothing method at the largest distances. Due 
to cancellation in the expression 



V B + Pf + x 2 ' 1 



dx 2 

this function would drop too fast. The two correlation lengths associated with T> B and T>f 
by the fit are almost degenerate and do not grow with temperature. 
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For the E correlators, a correlated fit of V E and Vf starts to be increasingly inadequate 
at temperatures T > T dcc . An indication is the E correlator amplitude A for T = 1.15 T dcc 
where the signal has an uncertainty of the same order of magnitude. Taking the parameters 
of this fit seriously, the most dramatic feature is the drop of the correlation length of V E 
across the phase transition. At (3 = 1.4 (T = 0.71 T dcc ), safely below the deconfining 
transition, we observe a weak splitting of the two correlation lengths, with A ~ 0.22 fm. 
At T = 0.9 T dcc we find still A ~ 0.18 fm. Immediately above the deconfining transition 
(at T = 1.02 T dec ) it has dropped already to 0.10 fm. 

B. Gauss-Exponential Fit 

For the purpose of the fit shown in Fig. 1, we assumed all four functions to have a form 
like A exp (-B \x\ — C \x\ 2 ) and required the same limit for \x\ — > for V E and Vf, as 
well as for V B and V B . 

In addition to this, have we considered remaining perturbative contributions to these four 
functions being proportional to l/|£c| 4 . For all temperatures and all correlation functions 
we found this contribution to be smaller than the signal by one or two orders of magnitude 
over the range of distances from 1 to 6 lattice spacings. Of course, also here only correlator 
data at distances d > 2 a have been included into the fit. In the following we therefor only 
discuss parameters which were determined by fits neglecting such a power-like perturbative 
contribution. 

The parameters A, B, and C are listed in Table III in lattice units. In the confinement 
phase, the coefficient C is relatively important only for T> E and T> B . While it remains roughly 
unchanged also at higher temperature for the magnetic correlator, it rises dramatically in 
the electric case in order to account for the rapid decay. On the basis of this fit, we have 
determined the structure functions T> E , Vf, T> B and T> B which are shown in Fig. 2 for all 
j3 values. 

In a first step the V 1 ' have been numerically obtained by integrating the difference 
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T> E,B — V E,B (as represented by the fits) from to |a;|. The integration constants V E ' B (0) 
are fixed by the condition that T>\ — > for \x\ — > oo. We find the functions V E,B positive. 
Finally, the T> E,B have been obtained subtracting T> E ,B from V E,B . 

The amplitudes V E (0), T>f (0), V B (0), V B (0) (which are related to the gluon con- 
densate at high temperature) and the integrated correlation lengths of T> E and T> B , 
ifnF = Io° d\x\V E ' B (x)/V E,B (0), are given in Table IV for the six temperatures. As long 
as T < Tdec the functions T> E,B are smaller by roughly a factor four than the corresponding 
T> E ' B . That means that self dual or antiselfdual fields are less dominant in this range of 



temperatures compared to zero temperature. The ratio T>i/T> has been found [22] consider- 
ably smaller below the deconfinement transition in an analysis of SU(3) correlators obtained 
with cooling. Since there are no cooling data for SU (2) so far in our temperature range, we 
cannot decide whether this is due to the other gauge group or whether there is a systematic 
effect of cooling to enhance the contribution of locally (anti) self dual configurations to the 
correlator at higher temperature but still in the confined phase. 

The condensates V E (0) and V B (0) only slightly grow across the transition and start 
to rapidly grow not before T > 1.2. ..1.5 Td cc - Up to these temperatures the T> E,B (0) have 
increased to a size comparable with T> E,B (0). 

The integrated correlation lengths of V E,B shrink in a rather continuous way over the 
temperature range considered from 0.3 fm (electric and magnetic) at T = 0.7 T dec to 0.13 
fm (magnetic) and 0.07 fm (electric) at 1.9 T^ ec (Fig. 3). The smooth change of the four 
structure functions with increasing temperature is shown in Fig. 2. 



V. DISCUSSION 

At first, we find our data approximately 0(4) rotational invariant also at T ^ as long 
we consider the confinement phase. The small difference between the corresponding electric 
and magnetic correlators (both for T>n and T>±) probably can be attributed to overlapping 
fields of (anti)calorons, within a small density expansion. This picture, however, breaks 
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down in the deconfined phase. In the confinement phase, the correlator V E practically 
coincides with V B (see Fig. 1 and note there the upward shift by a factor two in the latter). 
At the transition, the breakdown of D E (in the correlation length) is the most notable result 
of our analysis. 

Comparing the two methods of fitting we find that the condensates estimated from the 
constrained, purely exponential fits in terms of two structure functions V E and Vf tend 
to be bigger than the condensates from the second fit. This is because the latter better 
accounts for the real shape of the correlators near zero distance. 

The exponential fit also does not properly describe the correlator V B at large distances. 
This observation holds for all temperatures. Within the deconfined phase, the E correlators 
cannot be described at all, by an Ansatz in terms of exponential T> E and T>f functions. 
Thus, the rapid change of the exponential correlation lengths (in particular, of the electric 
correlation length) at the transition merely reflects the difficulty to reach a constrained fit 
within these assumptions. 

The Gauss-exponential fit proves to be much more robust over the whole range of tem- 
peratures. Without any further assumptions concerning the shape of V E,B and V E,B , such 
a fit gives access to these functions which indeed seem to be poorly described by exponen- 
tials. We have to describe them by integrated correlation lengths which cannot be directly 
compared with the exponential correlation lengths of the supposedly exponential fits. In 
general, the integrated correlation lengths turn out to be bigger than the exponential ones. 

In contrast to the exponential correlation lengths of the constrained exponential fits, 
the integrated correlation lengths associated to the correlation functions V E < B and V E ' B , as 
reconstructed from the second fit, change rather smoothly. As the result of these smooth 
changes with temperature, the correlation function T> E drops (in correlation length) at the 
deconfinement transition. 

The reconstruction of the absolute size of V and V\ structure functions (condensates) 
from the fits of T>n and T>^ becomes increasingly difficult at higher and higher temperature. 
Below the deconfinement temperature and even near to Td cc as a characteristic feature of 
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the data, it is well-recognized, that the D\ part of the electric and the magnetic correlator 
is considerably smaller than the T> part. 

The dominance of (anti) self dual fields, however, associated with this observation is al- 
ready weaker in our present results, at temperatures below the transition, than it has been 
found for zero temperature. It is also weaker than what has been found analyzing cooling 
data for quenched SU(3) near the deconfinement temperature. In Table V this statement 
can be checked for the so-called non-Abelianicity of the magnetic correlator. It is roughly 
80 % at T = 0.7 Tdec and drops to than 70 % at Td ec . In the deconfined phase k b amounts 
to approximately 50 %, with an increasing uncertainty. 

The positivity of magnetic and electric T>\ is a clear feature of the data. This applies to 
the correlators obtained in this work by smoothing as well as to the measurements using the 
cooling method (for quenched SU(3)). Within a caloron analysis (which successfully works 
only below Td ec for the S77(3) cooling data) there is no compelling reason to have T>\ < 0. 
This is unlike the instanton liquid description for T = 0. A caloron analysis of the present 
data will be given elsewhere |KJ . 

The sum of T>(0) and Pi(0) is given by the extrapolation to zero distance of both Py 
and T>± for the electric and magnetic correlators. This quantity is known without big 
uncertainty and is also presented in Table V in physical units for the six temperatures. The 
behavior of condensates and correlation lengths is sketched graphically across and beyond 
the deconfinement transition in Fig. 3. Quantitatively, the correlation lengths are not very 
different from the results of the cooling method, applied on 577(3) configurations. 

The present work dealt with SU (2) pure gauge theory instead of the phenomenologically 
more interesting case of S77(3). However it has shown an interesting pattern of temperature 
dependence of the magnetic and electric correlation lengths related to the different structure 
functions. While the basic structure functions T> and V>\ change smoothly across the phase 
transition, the most dramatic effect, the breakdown of results from the interplay. A 
simple picture according to which T> would turn to zero in the deconfined phase is not 
corroborated by our data. 
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It seems to be worthwhile to investigate the field strength correlators for SU(2) using 
the cooling method, too, in order to compare it with the smoothing method applied here. 
It would be even more interesting to implement the RG smoothing method for the gauge 
group SU(3). This would then allow for a direct comparison with the pioneering results 
already available thanks to the Pisa group. 
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TABLES 



w(i,j) 


J' = l 


J = 2 


J = 3 


j = 4 


i = 1 (plaquettes) 


1.115504 


-0.5424815 


0.1845878 


-0.01197482 


i = 2 (6-link loops) 


-0.01443798 


0.1386238 


-0.07551325 


0.01579434 



TABLE I. Weight coefficients of the simplified fixed-point action for SU{2). 



Type 


P 


T/Td cc 




Ao X Tdec 




Ai x T dec 


X 2 /d-o.f. 


E 


1.40 


0.71 


3.92(2) 


0.326(5) 


1.59(1) 


0.354(1) 


15.92 


E 


1.50 


0.90 


6.94(5) 


0.274(6) 


3.09(3) 


0.314(1) 


10.11 


E 


1.55 


1.02 


17.4(8) 


0.156(2) 


5.016(6) 


0.318(1) 


10.84 


E 


1.60 


1.15 


270(290) 


0.055(6) 


12.80(4) 


0.2520(3) 


27.33 


B 


1.40 


0.71 


3.50(1) 


0.3455(5) 


1.45(1) 


0.368(1) 


7.271 


B 


1.50 


0.90 


6.59(3) 


0.2920(6) 


3.04(3) 


0.310(1) 


4.526 


B 


1.55 


1.02 


7.03(6) 


0.2832(7) 


4.45(3) 


0.290(1) 


0.823 


B 


1.60 


1.15 


7.41(4) 


0.2760(8) 


6.51(4) 


0.273(1) 


0.294 


B 


1.70 


1.47 


12.1(1) 


0.2386(9) 


14.9(1) 


0.2272(8) 


0.448 


B 


1.80 


1.88 


23.2(3) 


0.2004(8) 


35.2(3) 


0.1853(7) 


0.660 



TABLE II. Condensates and correlation lengths obtained from exponential fits of the form 
V E,B (x) = Aq exp (— |aj|/Ao) and T>f' B (x) = A\ exp (— |as|/Ai) for various temperatures T. The 
data was fitted for lattice distances d > 2a. Type E refers to V E and Vf , type B to V B and Vf . 
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0.130 


vf 


1.80 


0.095(1) 


0.6212(9) 


0.0291(2) 


0.130 




1.80 


0.013(2) 


0.5002(7) 


0.0258(2) 


0.043 


vf 


1.80 


0.013(4) 


0.539(2) 


0.1227(8) 


0.043 
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TABLE III. Parameters from the fit to the raw data of the correlators shown in Fig. 1 in 

E B E B 

lattice units. The parameters A, B and C correspond to a respective Ansatz for Py ' and V ± ' , 
V(x) = A exp (-B \x\ - C \x\ 2 ). 



a 


T/T dcc 


V E (0)/Ti c 


2>>)/T d 4 ec 


^f(0)/T d 4 ec 


2>f(0)/T d 4 ec 


£int X -^dec 




1.4 


0.71 


2.00(5) 


2.0(1) 


0.54(6) 


0.52(3) 


0.44(3) 


0.42(2) 


1.5 


0.90 


3.9(2) 


3.7(3) 


1.1(1) 


1.2(3) 


0.38(1) 


0.34(2) 


1.55 


1.02 


4.9(4) 


4.4(2) 


1.1(6) 


2.0(2) 


0.35(2) 


0.25(3) 


1.6 


1.15 


3.6(4) 


4.4(6) 


3.1(4) 


4.2(7) 


0.27(4) 


0.24(4) 


1.7 


1.47 


8.4(7) 


8.1(16) 


5.2(8) 


5.5(23) 


0.25(3) 


0.15(3) 


1.8 


1.88 


20.1(57) 


22.0(20) 


9.9(59) 


19.0(24) 


0.20(2) 


0.10(3) 



TABLE IV. Condensates and integrated correlation lengths of V E,B for the fit presented in 
Table III for various temperatures T. All quantities are expressed in comparison to T dcc . 
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p 


P £ (0)+Pf(0) 

V / 1 IV/ 

[GeV 4 l 

L J 


V B (0)+V B (0) 

V I 1 IV/ 

[GeV 4 l 

L J 




a B 
[GeV 2 l 

L J 


1.4 


0.020(1) 


0.020(1) 


0.79(3) 


0.10(1) 


1.5 


0.039(2) 


0.039(4) 


0.75(5) 


0.14(3) 


1.55 


0.047(8) 


0.051(3) 


0.69(15) 


0.14(3) 


1.6 


0.053(7) 


0.068(10) 


0.51(41) 


0.15(5) 


1.7 


0.107(12) 


0.107(30) 


0.60(23) 


0.15(10) 


1.8 


0.237(92) 


0.323(34) 


0.54(64) 


0.27(19) 



TABLE V. Some physically interesting quantities (for SU(2)) in physical units. The magnetic 
string tension a B is given by ( JlOj ) in terms of T> B . For definiteness, we took T& ec = 0.298 GeV for 
pure SU{2) gauge theory. 
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FIGURES 




123456 123456 

Distance r/a Distance r/a 

FIG. 1. Raw correlator data fitted according to a generic form A exp (— £?|x| — C|x| 2 ). The 
corresponding fit parameters are reported in Table III. Errors are determined by the jackknife 
method for a block size of 100. The dependence on the block size was found to be marginal. For 
the clarity of the plot the data and curves shown for T>f_ and T>f include a factor two shift upward. 
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Distance [fm] Distance [fm] 

FIG. 2. V 1 ■ and T> 1 ' reconstructed from the Gauss-exponential fit of Fig. 1 and Table III, 
here presented in physical units. 
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FIG. 3. Condensates and correlation lengths as listed in Table IV, here expressed in physical 
units and shown vs. (3. 
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